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Abstract 

The properties of magnetized plasmas are always investigated under the hypothesis that 
the relativistic inhomogeneities stemming from the fluid sources and from the geometry it- 
self are sufficiently small to allow for a perturbative description prior to photon decoupling. 
The latter assumption is hereby relaxed and pre-decoupling plasmas are described within 
a suitable expansion where the inhomogeneities are treated to a given order in the spatial 
gradients. It is argued that the (general relativistic) gradient expansion shares the same fea- 
tures of the drift approximation, customarily employed in the description of cold plasmas, 
so that the two schemes are physically complementary in the large-scale limit and for the 
low-frequency branch of the spectrum of plasma modes. The two-fluid description, as well 
as the magnetohydrodynamical reduction, are derived and studied in the presence of the 
spatial gradients of the geometry. Various solutions of the coupled system of evolution equa- 
tions in the anti-Newtonian regime and in the quasi-isotropic approximation are presented. 
The relation of this analysis to the so-called separate Universe paradigm is outlined. The 
evolution of the magnetized curvature perturbations in the nonlinear regime is addressed for 
the magnetized adiabatic mode in the plasma frame. 
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1 Motivations 



The analyses of the large-scale galaxy distribution [TJ [2], of the high-redshift type la su- 
pernovae [31 H] and of the Cosmic Microwave Background (CMB) observables [5j [6] seem 
to converge, these days, on a concordance model sometimes called ACDM scenario where 
A stands for the dark energy component and CDM accounts for the dark matter compo- 
nent. The ACDM scenario is just the compromise between the number of ascertainable 
parameters and the quality of the observational data. The quest for a concordance lore 
is also able to shed some light on the presence of large-scale magnetic fields in nearly all 
gravitationally bound systems we observe. Since we do see magnetic fields today over large 
distance scales, it seems natural to scrutinize their impact on the CMB observables. This 
is the motivation of a program aimed at bringing the unconventional studv of magnetized 
CMB anisotropies to the same standard of the more conventional adiabatiqj paradigm (see 
[3 El El EE] and references therein) . While different approaches to the problem are certainly 
available [HI [T51 [131 [EH H5] (see [32] for a more complete list of earlier references) , the path 
followed in [7J E] led to the calculation of the temperature and polarization anisotropies in- 
duced by the magnetized (adiabatic and entropic) initial conditions. The parameters of the 
magnetized background have been estimated (for the first time) in [91 [10] by using the TT 
and TE correlation^ measured by the WMAP collaboration. The obtained results^ show 
that large-scale (comoving) magnetic fields larger than 3.5 nG are excluded to 95 % C.L. 
and for magnetic spectral indices n-Q = 1.6°;^. These determinations have been conducted 
in the context of the minimal mACDM where m stands for magnetized. The addition of a 
fluctuating dark energy background pins down systematically larger values of the magnetic 
field parameters (see [TU] for further details). 

The results obtained so far assumed the simplest setup for the inclusion of large-scale 
magnetic fields in the pre-decoupling plasma and it is therefore mandatory to scrutinize if 
the main assumptions of the analysis are consistently posited. For instance in [7J El CEO] (as 
well as in nearly all other approaches) it was assumed that magnetic fields do not contribute 
to the electron-photon scattering. In [TT] this assumption has been relaxed by explicitly 
including the magnetic field contribution in the electron-photon scattering matrix. 

In Refs. [3 El El HO], as well as in other perturbative approaches to the physics of 
magnetized CMB anisotropies, a common hypothesis is that the intensity of the magnetized 
background is sufficiently small to describe magnetic fields within the standard Einstein- 

3 The wording "adiabatic paradigm" refers here to the situation where the unique source of inhomogeneity 
prior to photon decoupling is localized in the standard adiabatic mode of curvature perturbations. This is 
the situation contemplated by the minimal version of the ACDM scenario. 

4 Following the standard shorthand terminology the TT correlations denote the temperature autocorre- 
lations while the TE correlations denote the cross-correlation between the temperature and the E-mode 
polarization. 

5 In [2 [HI [9j [10] the magnetic power spectrum and the magnetic spectral index ne is defined with the 
conventions employed for the scalar modes of the geometry where the scale- invariant limit corresponds to 1. 
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Boltzmann hierarchy where the curvature perturbations as well as the density contrasts all 
remain in the perturbative regime. In similar terms, perturbation theory is extremely well 
justified for the treatment of baryon acoustic oscillations because of the absolute smallness 
of the relative temperature fluctuations. As far as fully inhomogeneous magnetic fields are 
concerned, the perturbative assumption is a direct consequence of the closure bound applied 
to the magnetic field intensity. The natural question, in this context, seems to be: can we go 
beyond the standard perturbative expansion and scrutinize its properties in a broader and 
firmer scheme? Can we understand which are the possible nonlinear effects in a gravitating 
and magnetized plasma to a given order in the spatial gradients? One of the standard tools 
to improve (and partially resum) the perturbative description of the relativistic fluctuations 
of the geometry is the gradient expansion where the guiding criterion is not the absolute 
smallness of the given perturbation relative to its background value but rather the number 
of gradients carried by each term and defining the order of the expansion. The extension of 
the general relativistic gradient expansion to the case of magnetized plasmas immediately 
suggests a physical connection with analog approximation schemes adopted in the discussion 
of weakly coupled plasmas in flat space-time. 

The modest aim of the present analysis will then be to combine the general relativistic 
gradient expansion with the drift approximation (customarily applied in plasma physics) 
and to derive a system of equations where both expansions can be implemented in a unified 
manner. We remind that a known tool for computing the properties of cold plasma involves 
an expansion in the gradients of the magnetic and of the electric fields. The details of the 
scheme depend upon the problem at hand and this richness is reflected in slightly different 
terminologies such as the drift approximation [18] (mainly adopted in kinetic theory) and 
the guiding center approximation [TjJJ [2D] (often used in the physics of cold plasmas and fully 
ionized gase^J). 

Similarly, the general relativistic gradient expansion has been discussed in several related 
contexts since the papers of Lifshitz, Khalatnikov and Belinskii (see, e.g. [22] |23| |2"H [25] ) 
but never in combination with either the drift approximation or the guiding centre approxi- 
mation. More recently various applications of the gradient expansion to inflationary models 
have been studied by Tomita (271 [23] as well as by Deruelle and collaborators [29] EUJ EI]- 
The neglect of the spatial curvature and of the spatial gradients to zeroth order implies that 
the obtainable solutions are, in a sense, opposite to the ones customarily discussed in the 
Newtonian regime and this is the reason why they are sometimes named anti-Newtonian. 
The latter solutions are also relevant for the so-called separate Universe picture stipulat- 
ing that any portion of the Universe larger than the Hubble radius th but smaller than the 
physical wavelength on the perturbation will look like a separate unperturbed Universe. The 
gradient expansion has been also applied to the problem of the evolution of non-Gaussianities 
(see, e.g. [321 E3] and references therein) and to the problem of dark energy [3J] with the 

6 The expansion in spatial gradients is also relevant in nonlinear magnetohydrodynamics (MHD in what 
follows) in flat space-time |21) . 
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purpose of demonstrating that the spatial gradients cannot be responsible of present-day ac- 
celaration (see e.g. [35J and references therein). An instructive approach to nonlinear power 
spectra in Einstein gravity has been developed by Noh, Hwang and collaborators in a series 
of interesting papers [361 E3 EH] • Finally the gradient expansion can be also employed in the 
investigation of the so-called sudden (or quiescent) singularities which arise in the context 
of dark energy models with baotropic index smaller than the one of a cosmological constant 
(i.e. the supernegative equations of state) [39|. 

Even if some of the considerations developed in this paper will be rather general, it is 
useful to bear in mind the essentials of the pre-decoupling plasma which contains different 
components interacting both gravit at tonally and electromagnetically. To fix the notations, 
it is convenient to separate the total action of the system, i.e. S to t in three distinct parts 

>Stot 'S'grav -|- S CTn -\- "S'plasma; (I'-Q 

where S'grav and S em denote, respectively, the gravitational and the electromagnetic contri- 
butions 



S'grav + S em — I d X \J — g 



2£l 16tt afi 3v 



(1.2) 

,p 1D7T J 

note that g = det(g fiu ) and F a/3 is the Maxwell field strength; S^iasma contains all the dif- 
ferent components of the plasma which can be written, in the context of the vanilla ACDM 
paradigm, as 

Splasma = Sei + Sj + S u + Scdm + S\] (1-3) 

S ei denotes the contribution of electrons and ions; S 1 and S u are, respectively, the contribu- 
tions of the photons and of the (massless) neutrinos; S c d m and S\ account for the cold dark 
matter and for the dark energy. 

Before photon decoupling the condition of validity of the general relativistic gradient 
expansion and of the guiding center approximation are both verified and it is therefore 
extremely interesting to derive and study the evolution equations describing the plasma 
without assuming that the geometry is, a priori, conformally flat. Various questions can 
be addressed in this scheme such as, for instance, the corrections induced by the gradients 
of the geometry on the two-fluid plasma description, on the MHD reduction as well as on 
all the plasma processes which are relevant for a correct description of magnetized CMB 
anisotropies. Needless to say that the present results are also applicable in the case of other 
magnetized systems in curved space-time not necessarily connected to CMB physics. 

The present paper is organized as follows. In section [2] the decomposition of the in- 
homogeneous geometry will be introduced with special attention to the themes which are 
mostly relevant for the present discussion. In section [3] the evolution equations of weakly 
coupled plasmas will be derived in the case where the metric is not assumed conformally 
flat and, in particular, without resorting to the standard separation between background 
geometry and relativistic fluctuations. In section H] we will scrutinize the way gradients must 
be counted in a magnetized plasma when the background geometry is fully inhomogeneous. 
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Section discusses the anti- Newtonian solutions in the drift approximation. In section [B] 
the quasi-isotropic MHD solutions are presented. Section [7] contains the derivation of the 
nonlinear magnetized adiabatic mode whose linearized counterpart is one of the elements 
of the simplest version of the mACDM scenario. Finally section [8] contains our concluding 
remarks. In the appendix (divided in two parts, i.e. appendix |X] and |B]) the complementary 
results have been collected with the purpose of keeping the paper self-contained. 



2 Fully inhomogeneous geometry 

The standard treatments of magnetized CMB anisotropies as well as the discussion of mag- 
netized plasmas in the early Universe customarily assume that the geometry is separated 
into a homogeneous background supplemented by its relativistic fluctuations. To go be- 
yond the latter description, the fully inhomogeneous geometry shall then be described in 
terms of N, N l and 7^ denoting, respectively, the lapse function, the shift vector and the 
three-dimensional metric tensor, i.e0 



#00 
9°° 
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9ij = 
N*Ni 
N 2 



'liji 



-r 



-N,. 



9oi 



(2.1) 



The decomposition of Eq. (12.11) is well known and it is sometimes referred to as the ADM 
decomposition from Arnowitt, Deser and Misner who applied it to the Hamiltonian formu- 
lation of the general relativistic dynamics (see, e.g., (101 SI])- In the ADM variables the 
extrinsic curvature fQj and the spatial components of the Ricci tensor become: 



2N 



-drla + (3) ViiV,- + (3) VjNt 



Q (3)pm a(3)pm ,(3) -nm (3)pn 

171 ij 3 i m ij m 



(3) pm (3)pn 
jn in, 



(2.2) 
(2.3) 



where ^Vj is the covariant derivative defined with respect to the metric 7^', d T denotes 
a derivation with respect to the time coordinate r and ^TJ- are the Christoffel symbols 
computed from 7^. Note that T™j =^ Y™: but only in the case Ni = (see appendix [Al for 
details). It is now useful to introduce the basic logic of the gradient expansion by looking at 
the whole system of Einstein equations with generic electromagnetic and fluid sourceqj. In 
their contracted form the Einstein equations read 



T 

T v 8 V 

A* 2 A« 



T 



(2.4) 



7 The Greek indices will run over the four space-time dimensions while the Latin indices will denote the 
spatial indices. 

8 As mentioned in sectionQ]the species present in the plasma prior to photon decoupling can be interacting 
directly with the electromagnetic fields and this will demand a more detailed description of the sources as it 
will be discussed in section [3] 
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where R u is the four-dimensional Ricci tensor while and T£ denote, respectively, the fluid 
energy-momentum tensor (accounting for the the global contribution of all the species of Eq. 
(11. 3p ) and the electromagnetic energy-momentum tensor: 



T, 



(p + p)u^u u - pg^, 



T v = — 
" 4tt 



5 

4 



(2.5) 



where g^ u u^u u = 1 (see also Eq. ( IA.12j) of appendix|A]for an explicit expression of the latter 
condition in the ADM metric ( 12. ip ). If not otherwise stated the covariant derivatives with 
spatial indices listed below will always refer to the three-dimensional metric. In other words, 
to simplify the notation, ®Vj = Vj. Details on the inhomogeneous geometry of Eq. (12. ip 
can be found in appendix |X] and will be quoted whenever needed. Bearing in mind these 
conventions, the various components of Eq. (12 .4p are given by: 



d T K - NTrK 



ViK - V k K" 



V 2 N 



N£ 2 P 



Nil 



u 



d T K 3 - NKK] - Nr\ + VWN 



P~ P 



5f-{p + p)uiV? + % J 



(2.6) 
(2.7) 
(2.8) 



where, for sake of simplicity, the shorthand notation u 2 = ^UiUj has been adopted. Having 
chosen Ni = in the general equations of appendix El the geometry appearing in Eqs. (12. 6ft - 
( 12. 8 p is described in terms of 7 independent functions (i.e. 7^ and N). It is useful to keep the 
lapse function arbitrary for practical purposes such as, for instance, the matching with the 
perturbative treatment of the magnetized fluctuations in the conformally Newtonian gauge 
(see, e.g., appendix [B]). The momentum constraint (i.e. Eq. ( 12. 7(1 ) can be used, in the 
framework of the general relativistic gradient expansion, to obtain the velocity field in terms 
of the extrinsic curvature evaluated to the preceding order in the expansion. By combining 
the trace of Eq. (12.81) with Eq. (12.61) . the standard form of the Hamiltonian constraint can 
be readily obtained 



24 



p + V + (p + p) 



K 2 - TrK 2 + r 
The traceless part of Eq. (12. 8p can instead be written as: 

V 2 N 



u 



d T K] - NKK] 



Nr 3 , + V t V 3 N 



-5 



i 2 P N 



(P + P) 
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UiU J 



Si 



r 1 



(2.9) 



(2.10) 



where the barred quantities define the traceless part of the corresponding variable, i.e. 



K 

K 3 - —8 J 

1 3 



T 

T"? — Fp. 

1 3 l ' 



T 

V - —5 3 . 

1 3 1 



(2.11) 



Recalling Eq. f[Ql) and using Eq. fl24H it follows that TyK 2 = K\ K) + K 2 /3. Bearing 
in mind the explicit form of Eqs. (I2.6p - (l2.8p . to zeroth order in the spatial gradients, the 
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peculiar velocities as well as the spatial curvature are neglected. From the momentum 
constraint, the zeroth-order results determine the peculiar velocity which can be used as 
an input for the following order in the expansion. By iteration the first-order correction to 
the geometry can be determined. While the zeroth order of the linearized approximation 
is by definition homogeneous, the zeroth-order of the gradient expansion can well be fully 
inhomogeneous but does not contain any spatial gradient. Conversely, the first-order depends 
upon the spatial gradients and upon the spatial curvature and so on and so forth. 

If the gravitating system is a magnetized plasma (such as the one present prior to pho- 
ton decoupling) the inhomogeneities of the electromagnetic fields will necessarily modify 
the trajectories of the charged species. In this case the nature of the physical system com- 
bines inextricably electromagnetic and gravitational gradients. The tenets of the (general 
relativistic) gradient expansion must be consistently combined with the expansion in spa- 
tial derivatives usually adopted in the analysis of weakly coupled plasmas in flat space-time 
[T9| [20] both in the two- and one-fluid approximations. There are different ways of introduc- 
ing the guiding centre approximation and the simplest one is to think of a gradient expansion 
of the large-scale magnetic field, i.e. denoting with B the (flat-space) magnetic field we can 
write that 

B t (S, r) ~ Bi(S , r) + (x j - x^djB, + ... (2.12) 

where the ellipses stand for the higher orders in the gradients leading, both, to curvature 
and drift corrections. A similar expansion can also be written in the case of the electric field 
with the caveat that, in a plasma, electric fields are anyway screened for typical length-scale 
larger than the Debye radius: 

Ei(x, r) ~ Ei(x , t) + (x j - 4)djEi + ... (2.13) 

In a flat-space plasma to zeroth-order in the expansion of Eqs. (I2.12p and (I2.13P only the 
time derivative of the magnetic fields are kept. To first-order the spatial derivatives of the 
magnetic field can then be used as an input to deduce the electric fields. The first deriva- 
tives of the electric fields (obtained to first-order) can be used to deduce the second spatial 
derivatives of the magnetic fields and so on. If we ought to combine the general relativistic 
gradient expansion with the drift approximation, the essential step is the generalization of 
the two-fluid description and of the usual MHD reduction to the situation where the fully 
inhomogeneous geometry is parametrized as in Eq. (12. ip . This analysis will lead automati- 
cally to the correct dynamical variables whose explicit form can be compared, for instance, 
with the corresponding variables deduced in the context of conformally flat geometries which 
are usually assumed in more conventional perturbative expansions. 
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3 Charged fluids and gradient expansion 



The electromagnetic interaction affects the evolution of the electron-photon-ion system while 
it affects only indirectly the evolution of the weakly interacting species. Since large-scale 
magnetic fields gravitate, the relativistic fluctuations of the geometry are modified by their 
presence via Eqs. ( 12. 6p . (I2.7P and (12. 8p . The purpose of the present section is to deduce the 
usual magnetohydrodynamical (MHD) reduction to lowest order in the gradient expansion, 
i.e. by allowing an underlying geometry with potentially large inhomogeneities over suffi- 
ciently large length-scales. The obtained equations generalize previous perturbative results 
(see, in particular, |9]) and will be cross-checked in the appropriate limits. The Maxwell 
equations in a four-dimensional curved space-time can be written as 

V^ = 47rf, V^ = 0, (3.1) 

where V M is the covariant derivative defined with respect to the four- dimensional metric g^ v 
while F^ u denotes the dual field strength. In terms of the ADM decomposition of Eq. (12. ip 
the field strengths and their duals are: 

F° l = -jp, F^ = - lmkV ^—, (3.2) 

f2i Cm 

F m = -jp, F^ = lmk ^ k —. (3.3) 
The totally antisymmetric Levi-Civita tensor rfi* is defined as 

e ijk 

V ljh = Vijk = V^djki (3.4) 

v7 



where e^k is the Levi-Civita symbol in flat space. The Maxwell field strengths with doubly 
covariant indices are given in appendix |A] (see Eq. ( IA.27P ) and simplify when the shift vector 



vanishes. Defining with e the electron charge, the total current appearing in Eq. (13. ip is the 
sum of the currents of the electrons and of the iong^l 

J =J( C )+J(i), J(e) = -en e M e , j (i) =eniWi, (3.5) 

where n e and fii denote, respectively, the concentrations of electrons and ions. The generic 
four-velocity w M satisfies g^ v u^,u v = 1 which means that 

coshy i v l coshy 1 

u — , u = , cosh y — — ; =, (3.0) 

N N J 1 _ |L v ; 



9 In Eq. p.5j) as well as in other equations the (lowercase) roman subscripts label the species (for instance 
electrons or ions) while the (lowercase) italic subscripts (or superscripts) must be understood as (spatial) 
tensor indices. 
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where v 2 = ^ijV l vK Note that, from now on, we shall deal, without loss of generality with 
the case Ni = 0. For immediate convenience the notations of Eq. (13. 6p can also be recast in 
a slightly different form by defining v % = v l /v: 



u 



v % sinh?/, 



UiU J 



Vivi sinh 2 y, u l u = v l N cosh y sinh y, 



(3.7) 



where v l v^ij = 1. Bearing in mind Eqs. ( I3.5p - (l3.7p . Eq. ( 13. ip demands the validity of the 
following system of equations: 



di 
0, 



V7 r 



N 
N 



47re[ni - n e ], 



V T ni 

N 
0. 



jki gr, 



'fkr. 



Aire 



mv {i) - n e v (e) 



0. 



sfl pi 



(3.8) 
(3.9) 
(3.10) 



where the electron and ion concentrations have been rescaled as: 

«(i) =n- x ^Jl cosh?/;, n(e) = n e v/7cosh?/e- (3.11) 

From the covariant conservation of the electron and ion currents (i.e. V^j^ = and 
^nJa) = 0)) the evolution equations for the rescaled electron and ion concentrations 

d T n e + d k [n e v k e ] = 0, d T m + d k [n^\ = 0, (3.12) 



imply that, to lowest order in the gradient expansion, n e and n ; are constant. Equations 
( 13. 8 p , ( 13.91) and ( 13.101) assume a simpler form by an appropriate rescaling of the electric and 
magnetic fields: 



d ■ E = Ane[n\ — n c ], 
d T B + 8 x E = 0, 

rii Vi - n e v e 



d-B = 0, 



d x B = 4ire n, V\ — n P v P + d T E, 
where the rescaled electric and magnetic fields are given by: 

E i = Vl £ i B > 



N 



v 2 

N 



(3.13) 
(3.14) 

(3.15) 



(3.16) 



The shorthand notation employed in Eqs. ( 13. 13f) . (13. 14ft and ( 13 .151) implies for a generic 
vector A\ 



8 ■ A ee d l A\ (d x AY = dj 



N~/ ik y jn r) nmk A r ' 



(3.17) 



The vectors appearing in Eqs. (13. 13j) . (13. 14j) and ( 13 . 1 5|) become three-dimensional Cartesian 
vectors in the limit when the background geometry is homogeneous, isotropic and confor- 
mally flat, i.e. 



N -)> a(r), = a 2 {r)8i 



(3.18) 
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Using Eq. ( 13.1 8p into Eq. ( 13 .lip we have that, for ?/ c < 1 and Hi< 1, the electron and 
ion concentration^ become n e = a 3 n e and n\ = a 3 hi. Using Eq. (13. 18[) into Eqs. ( 13. 13ft - 
( I3.15P the standard differential operators are recovered, i.e. d — > V, d x A — > V x A. In 
general terms, the system is not conformally invariant as it can be argued by looking at 
the form of Eq. (I3.17P and by appreciating that indices are raised and lowered in terms 
of 7y. Various discussions obtained in the limit defined by Eq. ( 13. 18j) |JJ [TO] (see also 
[3 |8]) can be generalized to the fully inhomogeneous situation. Without dwelling on all 
possible generalizations we shall focus our attention only on a consistent implementation of 
the two-fluid and one-fluid descriptions. 

Neglecting, for a moment, the terms responsible for the momentum exchange between 
electrons and ions the covariant conservation of the energy momentum tensor of the charged 
species implies 

V M IJJ = Jl e) i^, V,T^ = (3.19) 

where 

Tfe) = Pe«, Tffi = (3.20) 

Using Eq. (lA.20j) and choosing the free index of Eq. (13. 19j) to be time-like, the explicit 
evolution of the energy density of the electrons can be obtained and it is: 



d T [p c cosh 2 ?/ e ] + N 2 d k 



cosh y e sinh y e v e 



-NKp e cosh y c + 



3V k N + NT- 



-NK, 



kj 



PeV 



k vl sinh 2 y e + U kj 



p e cosh y e sinh y c v k 

j { : ] £ k , (3.2i) 



where the contribution of the anisotropic stress of the electrons H l J has been also included 
for completeness. The same expression holds also in the case of the ions by flipping the sign 
of the electric charge and by replacing e — > i in the various subscripts (i.e. p c — > p;, y c — > y\, 
... and so on and so forth). By choosing the free index of Eq. ( 13 . 1 9[) to be space-like the 
resulting equation is 

d T [p c cosh y e sinh y e v k e ] + d m [p e sinh 2 y e v™v k e + U^ k ] 

—2NK m p e sinh y c cosh yevT + V fc A^p e cosh 2 y e — NKp e sinh y e cosh y e v k 

+ [V m N + NT n nm }[p c v™v k e ^ 2 y e + U km ] + NT k mn [pv m v n s 2 {y) + U mn ] 

= A e) ^ - lm^\ knm B^ + C ei + C e7 , (3.22) 

where C e i and C e7 denote the collision terms of the electrons with ions and photons. Similarly 
the evolution equation for the ion velocity field can be obtained by replacing e — > i in the 



10 Elcctrons and ions are non-relativistic across matter-radiation equality and around decoupling. Still, 
because of the masses of the electrons and ions, the conformal invariance of the whole system is broken (see, 



e.g. 01). 
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relevant subscripts in full analogy with what has been already suggested, after Eq. ( I3.2ip . for 
the evolution equations of the energy density. Equations (I3.2ip and (I3.22j) can be expanded 
in gradients and, to lowest order, the evolution of the electron and ion energy densities can 
be derived from Eq. (13. 21ft : 

d T p e — NKp c = —en e — — , d T p\ — NKp\ = en x — — , (3.23) 

7 7 

where, generically, the following notation will be employed throughout the paper: 

F ■ G = ~f mn F m G n , (F x G) k = lilMF n G m r) ak . (3.24) 

The reduction from the two-fluid to the one-fluid description follows the standard steps of 
the flat-space time case (see, e.g. [12]) but with the difference that the terms stemming 
from the fully inhomogeneous nature of the underlying geometry will be consistently taken 
into account. The idea of the one-fluid reduction is to pass from a description holding for 
the two (or more) separate species to a one-fluid system where the dynamical variables are 
global, like the total current, the baryon energy density, the center-of-mass velocity of the 
electron-ion system and so on and so forth. The sum of the ion and of the electron energy 
densities is usually referred to as the baryon density and its evolution equation is obtained 
by summing, term by term, the two equations appearing in Eq. ( 13.231) : 

d T p h - NK Ph = — , (3.25) 

7 

where pb and J are, respectively, the baryon density and the total current: 

Pb = Pe + Pu J = e^riiVi - n e v e ). (3.26) 

In Eq. (13.251) . the contribution of the electric field has been kept since it multiplies the total 
current which will turn out to be related with spatial gradients of the magnetic fields as well 
as with the spatial gradients of the geometry possibly multiplied by terms containing the 
magnetic field itself. From Eq. (I3.22|) the evolution equation for the electron velocity can 
be written as: 



d T v k e + Nd k N - Q k v{ = r: "" A " 



/W7 



E k + {v e x B) } 



NT^v k -v k e ) + - P -lNT cl {v k -v k ), (3.27) 

" Pe 



where 

-5 k + 2NK 



g; 



N 



(3.28) 
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In Eq. (I3.27P the collision terms have been included and the shorthand notation of Eq. 
(I3.24p has been used for the vector product. Similarly, the evolution equation for the ion 
velocity is given by 

d T v? + Nd k N-g*vl . n ' ]X2f 



(3.29) 



PW1 



E k + (vi x B) k 



+ NTUv k e -v k ) + ^NT h (v k -v k ). 



By summing up Eq. (13.271) (multiplied by the electron mass) and Eq. (I3.29P (multiplied by 
the ion mass) the evolution equation for center of mass velocity of the electron-ion system 



m e v k + m\v k 



m e + m\ 

can be obtained and the evolution equations for the baryon-lepton-photon system are 

4 



<9 r p 7 = 


4 

-KNp 


d T v k = 


qU- 


d T v k = 


k-- 



;Nd k 



N v vr 



(JxB) k N 2 4 p*y . k fcN 



ryp h (l + m e /mi) 3p b 

N \ ( mk 



Nd k N + NT,Jvt 



(3.30) 
(3.31) 

(3.32) 
(3.33) 



where v k and p 7 denote, respectively, the photon velocity and the photon energy density. 
The possibility of describing the baryon-lepton-photon fluid as a unique physical entity is 
a direct consequence of the fact that the electron-ion collision rate is much larger than the 
electron-photon (or ion-photon) rate. While the electron-photon rate increases with the 
temperature the Coulomb rate decreases. The meeting point of the two rates occurs close 
to the MeV. Equations (I3.3ip - fl3.33p can then be used below the meeting point of the two 
rates. 

Equation f!3.29j) (multiplied by en{) can be subtracted from Eq. (I3.27P (multiplied by 
en e ) leading to the generalized Ohm equation, i.e. the evolution equation of the total current: 

dJ k 



3r - e > J '— e 



v[ dj(niv() - v e dj(n e v e 



- e(ni - n e )Nd k N 



1 E — erii — + en e — 



Air 

+en e NT l , i [ I 
-f 



Pi 



Pe- 



rn, 



1lh 



(Wj -n e )(m e + mj) k _ (m t + m e ) k 
m;n e + n;m e e(nim e + n e mi) 



e 2 n e rii(m e + m-^N 2 
eN 2 



(nim e + n e rrii) \ 1 m\ e m, 



(J x By 
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i f / / ^ i 7 ^7 
+ Q e ^l 

3 I V mi m e 



m e (rriin e + m e n{) mi{m\n e + m e n\) 



e7 



e{m- x n & + nim e ) Vm e / e{m{a c + nim e ) 



+ 
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(3.34) 



where Eq. (I3.12p have been used and where the plasma frequencies for electrons and ions 
are defined, respectively, as 



4ne 2 n P N 2 



pe 



47re 2 niiV 2 



(3.35) 



m eA /7 p m iA /7 

Since the plasma is globally neutral, n e = n; = n (with0 n = rj-^n^). Thus Eq. (13.351) 
greatly simplifies also because of the smallness of the ratio between the electron and ion 
mass: 



dJ k 



+ 



NT ie 5* + 

J 6n m e ■ 



00 



pc 



47T 



E k + 



k N 
v h x B) H dj ( pe7 

en 



Vidj(nii4) - v^dj(n e v e 
(fx B) k - 



en Q 



(3.36) 



In Eq. (I3.36P the terms containing the spatial gradients have been kept to illustrate the 
analogies and the differences with the customary MHD discussions in flat space-time. The 
third and fourth terms appearing in the square bracket of the second line of Eq. ( I3.36P 
are, respectively, the thermoelectric term (containing the gradient of the pressure) and the 
Hall term (containing the vector product of the current and of the magnetic field). The 
thermoelectric term contains pressure gradients and it is therefore of higher order in the 
gradient expansion. The terms containing the velocity field multiplied by its own gradient 
are even smaller than thermoelectric term. The Hall term is of the same order of the 
thermoelectric term. Equation (13.361) can also be schematically written as 



d T J k + (NT5f-gf)J l = S k 



(3.37) 



where T = T ie + [4p 7 /(3n m e )]r e7 ) denotes the sum of the electron-ion and electron photon 
rates. Equations (I3.32p and (13.331) can be combined by noticing that, for sufficiently early 
times prior to decoupling, the baryon and photon velocity coincide. The resulting equation 
generalizes the standard evolution equation for the baryon-photon velocity which plays a role 
in the semi-analytic treatment of the magnetized CMB anisotropies at small angular scales 
[H] . Another possible generalization concerns the propagation of electromagnetic waves in the 
plasma which is relevant for the Faraday rotation of the polarization plane of the CMB (see, 



11 Recall that rjh denotes the ratio between the baryonic concentration and the photon concentration and it 
is given by ?7b = 6.219 x lCT 10 (/iQS!bo/0.02773)(2.725 K/T 7 o) 3 where ho is the indetermination on the present 
value of the Hubble rate and fibo parametrizes the present value of the critical fraction of baryons; T 7 o is 
the CMB temperature. 
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e.g. [44J). In conclusion all the evolution equations customarily employed for the description 
of weakly coupled plasmas in the linearized approximation can be generalized to the case of 
a fully inhomogeneous geometry expressed in ADM variables. 

The evolution equations derived in this section reproduce, in the conformally flat limit, 
the standard perturbative results of Ref. [9]. As an example consider Eqs. ( I3.3ip . ( I3.32p and 
( I3.33P describe the evolution of the photon-baryon system to lowest order in the gradient 
expansion. All the quantities appearing in the latter equations depend both on x and r 
and generalize the perturbative approach to the study of the relativistic fluctuations of the 
geometry in a cold plasma. To reproduce the perturbative results it suffices to identify: 

N(x,t) = a(r)[l + (p(x,r)], 7^(2, r) = a 2 (r)[l - 2ip(x, r))5 ih 

p 7 (f,r) = p 7 (r)[l + 5 7 (x,r)\. (3.38) 

Provided 0, ip and £ 7 are perturbatively well defined, the choice of Eq. ( I3.38|) corresponds 
to the case of the conformally Newtonian gauge discussed, for related reasons, in appendix 
E Using Eq. ( ET5B1) . Eqs. ( ET5Tj) . (13T32D and ( 13331) then become: 

<9 T 5 7 = 4d T ij - ■ Vry, (3.39) 
<9 T Uy = — ^V5 7 — V(f) + aTjeivb — Vry), (3.40) 

d T v h = -Uv h + J j — - V0 + ^ar 7e (u 7 - v h ), (3.41) 

which coincide with the system discussed in [9]. For instance Eqs. (13.391) . (I3.40p and (I3.4ip 
have been used to derive the appropriate initial conditions for magnetized CMB anisotropies 
in the tight-coupling approximation. In connection with Eq. (I3.3ip we should mention a 
general feature of the interplay between the gradient expansion and the more conventional 
perturbative expansion leading to Eqs. (I3.39p . (I3.40p and (13.41 p : the term containing the 
velocity field in Eq. ( 13 . 3 1 [) is of higher order in the gradient expansion but it is necessary 
to recover the correct linearized result. This occurrence simply signals that the gradient 
expansion and the standard perturbative expansion do not commute: the linearization of 
the equations written at the lowest order in the gradient expansion leads to a set of equa- 
tions which is different from the one obtained by linearizing the exact equations. A useful 
compromise, as shown in Eq. ( 13 . 3 1 j) . is to keep in the lowest order of the gradient expansion 
all those terms leading, after the linearization, to the standard perturbative equations in a 
given gauge. 
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4 Counting gradients in weakly coupled plasmas 



Let us consider, for sake of concreteness, typical length scales of the order of (or larger than) 
the Hubble radius (and hence larger than the Debye length) and typical time derivatives 
much smaller than the plasma frequency (see e.g. Eq. (I3.35P ). In weakly coupled plasmas 
the plasma parameteio g p i a sma is, by definition very small 

1 A , [cm , 7 f h 2 SL hQ \V2 

ft— = wl = ^V^^= 2 ' 308 x 10 ^(™) • 

v - = ad = < 4i) 

where the reference temperature T has been taken to coincide with the photon temperature 
(recall that the electron and proton temperatures coincide up to a small quantity which 
is the ratio of r-^ to the electron-photon rate [Tl [8]). Given the smallness of (y'piasma, 

the 

electron-ion mean free path is much larger than the Debye scale while the corresponding 
collision frequency is much smaller than the plasma frequency, i.e. 

r ei In Ac 

T— - — in A C , - ffplasma, l 4 ^J 

A ei 487T LU pe 24V27T 

where A c = (lSy^/fi'piasma) is the argument of the Coulomb logarithm. Since r H greatly 
exceeds A e i the space-time gradients of the system under consideration will be ordered as 
follows: 

V<0 U <<0 U <<0 U> ar<fKJ. (4.3) 

where V and d T denote, respectively, the magnitude of a typical spatial gradient and the 
magnitude of a typical time gradient. Let us then consider Eq. (14. 3D in conjunction with 
Eq. ( I3.37P and with Eq. (I3.15P whose explicit form, in terms of the total current, can also 
be written as: 

x B) k = 4irJ k + d T E k . (4.4) 

Because of the hierarchies established in Eq. (14.31) the time gradient of the total current 
must be negligible in comparison with the term containing the total rate. Similarly the total 
current must be much larger than the displacement current, i.e. 

d T J k ~ o(g k .A < NTJ k , AnJ k > d T E k . (4.5) 



12 The definition of the Hubble radius rn refers, strictly speaking, to the homogeneous and isotropic case 
and can be slightly ambiguous in the fully inhomogeneous situation which is the one treated in the present 
paper. For the moment, rjj must be considered as the inverse of the trace of the extrinsic curvature (possibly 
evaluated in a specific gauge). 

13 The plasma parameter quantifies, by definition, the inverse of the number of particles present in the 
Debye sphere, i.e. the sphere whose radius is given by the Debye length. 
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Equations (13.371) and (14. 3 p together with Eq. (14.51) also imply the following pair of inequali- 
ties: 



d T J k < u 2 pe d T E k , 



d 2 T J k < oj 2 pe J k . 



(4.6) 



But then Eq. (14. 6[) means that the Ohm equation of Eqs. (13.361) and (13.371) reduces to the 
following Ohm law: 



J 



a 



a 



pc 



(4.7) 



where the conductivity a is not bound to be homogeneous as in the case of previous treat- 
ments [9]. Equation (14. 7p generalizes to the fully inhomogeneous situation the standard 
result of the conformally flat limit already mentioned in Eq. (I3.18P (see, e.g. [8], third refer- 
ence). Note that, in the latter limit, the spatial derivatives of the extrinsic curvature and of 
the determinant of the metric are all vanishing; moreover the total current, the electric and 
magnetic fields, the conductivity are all rescaled through different powers of the scale factor 
with respect to their flat-space values. Consider now Eq. (14. 7p written in its explicit form. 
Because of Eqs. ( I4.4p . (14. 51) and (14.61) the Ohmic electric field can be expressed as: 

1 



E' 



AlXO J 



Nf k ^ n Vmnk B r ' 



N 



qpi v[ B 9 , 



(4.8) 



where the second equation reported in Eq. (13. 13f) has been used to simplify the obtained 
result. Equation (14. 8f) still implies, for weakly coupled plasmas, that in the drift approxi- 
mation the electric fields depend upon the gradients of the magnetic field in the baryon rest 
frame. This means that the electric fields are higher order in the gradients. From Eq. 
the generalized magnetic diffusivity equation can be derived and the result is 



47T 



N 



9t b* + —d a - 7 fes 7 r \„a 



a 



NY k ^ n Vnmk B r 



Or 



i hs i ru Vsua fq i 9P vlB° v ^ 



(4.9) 

where, for immediate convenience, we preferred to avoid the shorthand notations employed 
before. The analog equation for the electric field is instead 



d T E i 



-d T 



rjipi 



+ 



di 

4:71(7 J 



dA[N^ k ^ n nmnk W 



(4.10) 



In the conformally flat limit mentioned in Eq. (13 . 1 8j) . Eqs. ( 14. 9 p and (14. 10|) coincide, respec- 
tively, with their flat space counterparts, namely: 



dB 

— = V x (v h x B) 
or 



1 



Ana 



-V 2 B, 



dE 



~% x B) 



V 2 E, 



(4.11) 
(4.12) 



where, according to Eq. (13 . 16}) . B = a 2 B and E = a 2 £. 



16 



Let us now rewrite Eqs. (12. 6p . ( 12. 7p and (I2.8P in a more explicit form which will turn out 
to be useful in the forthcoming sections: 



d T K - NTrK 2 + V 2 N = Nij 



3p + p 



+ (p + p)s 2 (y) + pb + pe 



{p + p)v iS {y)c(y) - -^-r] mnp j im j jn j pk B k E j 



d T K{ - NKKj - Nr{ + V;V j iV = Ni 2 P 



p- p 



$1 ~ (p + p) s 2 (y)viV J 



(4.13) 
(4.14) 

(4.15) 



where the notations of appendix [S] have been used (see, in particular, Eqs. flA.25p -f fA.26P 
and f lA.32p -f lA.34p ). To zeroth order in the gradient expansion the terms containing s(y) 
are subleading since y = v/N and v 2 = jijV V. By combining Eq. (I4.13P and the trace of 
Eq. (I4.15P the terms d T K can be eliminated and the resulting expression is the Hamiltonian 
constraint 



K 2 - TrK 2 + r = £ 2 P [2(p + p)s 2 {y) + 2p+ (p E + 3p E ) + (p B + 3p B )] 



(4.16) 



coinciding, as expected, with the expression already obtained in Eq. (I2.9p in the light of the 
values of pb and pe- We recall that the velocity field, the energy density p, the pressure p 
are global quantities given by the sum over the individual species: 



(4.17) 



and obeying conservation equations which can be obtained, species by species, using the 
results of Eqs. (1A.20I) and (1A.21|) reported in appendix [At 



d T p + Nd k 



(P + P) v k 
N 



NK{p + p) = 0, 



dr[(p + pY\ + N 2 d k (pY k ) + N 2 d k W l + (p + p)Nd t N - 2NK)(p + p)v 

-nk(p + P y = 0. 



(4.18) 



(4.19) 



As already mentioned after Eq. (14. 8p . the baryon rest frame is particularly useful for the 
treatment of the finite conductivity effects. In a perfectly conducting medium (i.e. a(x, r) — > 
00) the Ohmic electric field is perfectly screened. Owing to this occurrence, it is customary 
to define the plasma frame where the electric fields are set to zero. In the usual perturbative 
expansion defined, for instance, in appendix (Bj the plasma frame coincides with the baryon 
rest frame. In the fully nonlinear case, however, the two concepts do not necessarily coincide. 
It is finally appropriate to mention that there are also nonlinearities associated with the Ohm 
law itself (see Eq. (I3.36P ) and an example along this direction is the so-called nonlinear Hall 
effect. This effect comes by retaining the terms J x B in Eq. (I3.36p . The Hall term leads, 
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in the magnetic diffusivity equation (14. lip , to a term of the type V x [B x (V x B)]/(en ). 
The nonlinear Hall term can partially balance or even become greater than the dynamo 
term, under certain conditions [43J . The present framework paves the way for the consistent 
treatment of the gravitating counterpart of the nonlinear effects typical of cold plasmas (see, 
for instance, [ZT]). 



5 Anti-Newtonian drift approximation 

Consider, for simplicity, the case where N = N(t)\ then jij(x, r) corresponds to 6 unknown 
functions and all the terms containing at least one spatial gradient of N vanish exactly. The 
spatial curvature is neglected since it contains two spatial gradients of 7^. Bearing in mind 
the definition of Eq. (12.111) . Eqs. (14.131) and (14.151) can be written as 

Nil 



d T K = NTtK + 



2 , iV< -p 



(3p + p) - 2p A + 2p B 



(5.1) 



3f 2 N 

d T K = NK 2 + - p - 2p A - 2p B ], (5.2) 

d T K\ = NKK\ + 4 N nf , (5.3) 
d T p = NK(p + p), drB 1 = 0. (5.4) 

Concerning Eqs. ( I5.ip - (j5.4p few comments are in order. The Ohmic electric fields have 
been neglected since they are of higher order and, for sake of simplicity, 11^ denotes the 
total anisotropic stress. Note that, in Eqs. (I5.ip - (l5.4p . the total pressure and the total 
energy density have been separated, respectively, as p = p — p\ and p = p + p\ where 
Pa parametrizes the dark energy density contribution; p and p are the pressure and energy 
density of an ordinary fluid characterized by a (possibly inhomogeneous) barotropic index 
w = p/p. The approximations leading to Eqs. (I5.ip -( |5~4|) define the anti- Newtonian limit. 
Equations (15. ip and ( 15. 2ft can also be written in more explicit terms by introducing the 
barotropic index and by recalling that, as in Eq. ( 12. lip . K\ = K{ — (K/3)5j: 

3iV£ 2 

d T K = NK 2 + —^-[(w - l)p - 2p A - 2p B ], (5.5) 

d T K = jK 2 + NTtK 2 + ^[(3w + 2p A + 2p B ]- (5.6) 

By eliminating p between Eqs. (15. 5p and ( 15. 6p and by defining dt = N(r)dr, the resulting 
equation is given by 

w + 1 2 3(^-1)^—2 34 u ^ (3w-l) . , r n . 

d t K - ~^-K 2 + - TiK = —^-[(w + l)p A + 3 > b]- (5.7) 

The explicit solution of Eq. ( 15. 3p . Kl(x,t) implies that 



^(f,t) = -^=+4-r==, <4frt)= / x -u.i')ii:(r.r) ( ir. { :,.s) 

yi(x,t \Jl{x,t) 
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By now introducing the rescaled variable Ai(x,t) = [y(x, i)]( m+1 " 4 , and by recalling that 
PB&t) =p B (x)7" 2/3 , Eq. dSZZD 

<9 2 .M = l(w 2 - l)M^[Tr\ 2 + t P Tra 2 + 2f P Tr\cr] 
8 

+ 3( W + 4 1)24 Pa(^){i + i^MZ)M-^^} M, (5.9) 

where Ab(x) = Pb{x)/p\{x). Defining the following auxiliary coefficients 

C(x,t) = l(w 2 - l)[TrA 2 + 4Tm 2 + 24TrAa], V{x) = 3 ^ + /^ p A (£), (5.10) 
8 4 

Eq. (I5.9P becomes 

d 2 .M = C(x,t).M^ + 3 3 ^ ~ * A B (x).M- 8 /' 3 ( m+1 )'} X. (5.11) 

The general form of Eq. ( 15. lip does not have explicit analytic solutions. In various limits 
it is possible to integrate it once with respect to t in terms of space-dependent integration 
constants. The result, however, cannot be further integrated (or inverted) except that in few 
special cases. Eq. (15. lip can certainly be integrated numerically but this study is beyond 
the scopes of the present discussion. 

It is useful to look at Eq. (15. lip by bearing in mind the usual assumptions of the separate 
Universe picture stipulating that any portion of the Universe that is larger than the Hubble 
radius rn but smaller than the physical wavelength on the perturbation will look like a 
separate unperturbed Universe. Such a framework is justified in the case K\ = 0. In the 
latter case, the extrinsic curvature can be really thought as the inhomogeneous generalization 
of the Hubble parameter. Assuming, for consistency with this hypothesis, that large-scale 
electromagnetic fields are absent, Eqs. (l5.ip -( l5T3"j) formally coincide with the Friedmann- 
Lemaitre equations. Consider next a slightly more complicated situation, namely the one 
where pa = but K\ ^ 0. In this case Eq. (15.1 ip can be explicitly integrated and M. can 
be determined. Suppose, for sake of simplicity, that w = and that A B (x) = 0. In this case 
Eq. (15. lip can be solved exactly: C(x,t) becomes independent of t since al(x,t) — > 0. The 
determinant of the metric becomes, in the cosmic time coordinate t: 



7(f, t) 







2 3 TrA 2 1 




t - ti(x) 






8 H 2 J 



(5.12) 



where Hi is a space-time constant while, as already mentioned, TrA 2 can have an arbitrary 
spatial dependence. From Eq. (I5.12p we can argue that even if K\ ^ 0, TrA 2 affects the final 
solution in such a way that it can be reabsorbed in the definition of the initial time coordinate 
or, put it in different terms, K\ ^ but goes quickly to zero if the Universe expands since 
K\ — Af / ^/j. If the magnetic anisotropic stress is taken into account the full system must 
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be solved consistently and anti-Newtonian solutions do not seem sufficient even if general 
solutions of Eq. (15. lip would be available. Therefore it seems useful to explore a slightly 
different strategy and solve the system in the quasi-isotropic limit where the contribution 
of the intrinsic curvature and of the magnetic anisotropic stress can be explicitly taken into 
account order by order. 



6 Quasi-isotropic MHD solution 



In the MHD description, the fields E\ B l and J 1 defined in Eqs. (I3.16P and ( I3.26j) are 
solenoidal, i.e. 

diE* = 0, d l B i = 0, diJ* = 0, (6.1) 

so that flat space MHD is recovered in the conformally flat limit of Eq. (I3.18p . Let us 
consider, as first example of the present section, the following parametrization of the spatial 
geometry 

7ij (x, r) = a 2 (r)e~ 2 * (x> % + 2h tJ (x, r)), h\ = 0, N(r) = a(r), (6.2) 

where hij accounts for terms containing more than one spatial gradient. Using Eq. (16.21) 
into Eq. (12. 2p the explicit form of the extrinsic curvature becomes 

'd T *_ _ W 
a a 



d h 3 



(6.3) 



The conformally flat limit of Eqs. f]6.2[) and ( 16.31) corresponds to ^ — > and hij — > 0; 
following the same logic the total fluid pressure and the total energy density can be separated, 
respectively, as 

p(x,r) =p {0 \r) +pW(x,r), p(x,r) =p (0) (r) + pW(f,r), (6.4) 

where p^(x, r) and p^'(x, t) vanish in the conformally flat limit. Inserting Eqs. (I6.2I) - (I6.3I) 
into Eqs. f|4.8D - (l4.9l) the Ohmic electric field and the magnetic diffusivity equation are, 
respectively 14 !. 

1 



E 



Ana 



V x {e*B) -e-*4 x B, 



d T B + V x 

Ana 



e*V x \e*B 



V x {y h x B). 



(6.5) 
(6.6) 



The terms containing have been neglected since they are of higher order. Inserting Eqs. 
(16. 2p and (16. 3p into Eq. (12.61) the following relation can be obtained 

p + 3p 



+ 



a 2 f P 



(p + p) 



6 



,4* 



a 2 e" 2 *" 2 



s\y) + 



24vra 4 



E 2 + E * + B 2 + B 2 



(6.7) 



1 The vectors appearing hereunder and in the remaining part of the present section are the standard 
three-dimensional vectors. Similarly the vector products are standard vector products. 
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where the following notations have been adopted: 



E 2 = Sij&E', 



E 2 = 2h ij E i E i , 



B 2 = 5 ij B i B j , 



B 2 = 2h ij B i B j 



u 



2 -2*r 2 

a e 



s 2 (y) + 2c 2 (y)h ij v t v J ] } y 



v 2 = SyvV. (6.8) 



Note that E 2 <C E 2 and B 2 <C B 2 since E and B contain and are therefore of higher 
order in the expansion. Furthermore, in the resistive MHD limit (i.e. a ^> 1) pe *C Pb since 



Pe 



8na 4 I 16tt 2 (J 2 

^4* 



|V x 5| 2 + |V# x B\ 2 + 2(Vtf xB)-(VxB) 



2na 



[v h xB)-(VxB) + (v h x B) ■ (W x S) 



+ e 2 *|4 x B| 



(6.9) 



Equation ( j6. 9f) shows, in practice, that at the nonlinear level the baryon rest frame and the 
plasma frame are two complementary concepts. While in the baryon rest frame the electric 
energy density vanishes as a~ 2 for large conductivity, if the conductivity is not large the 
contributions of the baryon velocity and of the gradients of the geometry are of higher order 
in the spatial gradients in comparison with the magnetic contributions. The resistive MHD 
limit implies, as usual, the largeness of the magnetic Reynolds number defined as 



Rr, 



u 



> 1, 



1 



V 



Ana ' 



(6.10) 



where r] is the magnetic diffusivity, V" 1 denotes the typical length-scale associated with 
spatial gradients and u is the modulus of the typical velocity of the plasma element. The 
velocity field itself is of higher order in comparison with p, p, p-Q since, from the momentum 
constraint of Eq. ( I4.14p we can deduce 



U°Uj 



2 3* ^ 

2did T m + d k d T h1 + =f^r(£ x B) 



Ana? 



(6.11) 



where, following the notations of Eqs. (13. 7ft . and to lowest order in the gradient expansion, 
the expression vPui can also be written as u°Ui = e~ 2 ^ViC 2 (y). The term E x B can be 
written by using Eq. f 16 . 5 j) : 



Ex B 



(V x B) x £? + (W x B) x B 
Ana 



e-*(v h xB)xB. 



(6.12) 



Neglecting the Ohmic electric fields, Eq. (I4.15P leads to the evolution equation for h\: 



d 2 r h\ + 2Ud T h\ 



2* 



+ <9^<9%-^(W) 2 
3 



ir 



a 4l ~(P + P) UiU J - —51 + 



,4* 



Ana 4 



BiB* (6.13) 
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where, using Eq. (I6.2p into Eq. (I2.3P and keeping only the leading contribution in the spatial 
gradients, r^x, r) becomes, in explicit terms, 



(6.14) 



where V 2 \l/ = S^didj^ is the flat-space Laplacian and, similarly, (V\I/) 2 = S^d^d^ . Fi- 
nally, the trace of Eq. (I4.15P leads, in the case of Eq. (16. 2p . to the following condition 



<9 2 ^ + md T ^> - (2H 2 + d T H) = -e 2 * 



2V 2 ^ - (VVf 



2 (p-p p + p 



4* 



U — 



B + B £ 



(6.15) 



2 3 " 247ra 4 

Using the parametrization of Eq. (16.41) . Eqs. ( 16. 7ft and ( 16. 15ft imply the following relations: 



<9 2 ^ + Hd r ^ = 
<9 2 ^ + 5^9 r ^ 



6 

2e 2 * 



(3p (1) + p (1) ) + 



a 2 4 



-PB, 



[2V 2 ^ - (V*) 2 ] + a 2 



3?/ 2 = 4a 2 p(°), 2(?/ 2 - 9 r -H) = 4a 2 (pW + p(°>). 



(6.16) 

(6.17) 
(6.18) 



Equations (16.131) . (16.161) and (16.171) can be directly solved and the result, in the plasma 
frame, reads: 



q(x) + 



2 a 3w+1 e 2q( - g) 



3(3w + 5)H? 

2 q3w+1 e 2q(x) 



V 2 g(f) 



4 a 3 *"- 1 e 4 ^ 
60n(3w + l)H( 



B 2 (x) 



(6.19) 



{3w + l)(3w + 5)F 2 
4 a 3 " 1 " 1 e 4 <^ 



di&>q{x) 81 



+ d iq {x)d j q{x) - ^^^ Sl 



2ir(3w - l)(3w + 1) 



Bi(x)B J (x) - 



B 2 (x) 



6! 



(6.20) 



where g = g(x) is a generic function of the spatial coordinates encoding the dependence of the 
large-scale inhomogeneities. Note also that we used if; = "Hi/a; at the initial reference time 
T{. Equations (16 . 191) and ( 16.201) show that the relative importance of the spatial gradients 
does depend upon the barotropic index. 

The results obtained in the particular case of Eq. (16.21) can also be deduced in a more 
general parametrization, namely: 



lik = a 2 (r)[aik(x) + f3 ik (x, r)), 7 



a\r) 



a kl {x)-p kl {x,r)l N{r) = a{r), (6.21) 



where (3ij(x, r) contains the contribution of the gradients while ctij(x) is fully inhomogeneous 
but does not contain any gradient. The conformally flat limit of Eq. ( I6.2ip does correspond 
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to a.ij — > 5ij and ^ — > 0. The indices of are raised and lowered by using a^. According 
to Eq. (12.21) . the extrinsic curvature and its contractions become: 



AV = -i(nS( + ^^j, K=-^(3H+^d T /3^, TrK 2 = L[3U J + U0 r > ). (0.22) 



where, as in the previous case, H. = <9 T lna. From the momentum constraint it also follows 
that 

V k d T tf - V t d r p = 2af P v i c(y)s(y)(p + p). (6.23) 

The explicit form of the momentum constraint suggests to look for the solution in a separable 
form, namely, 

0{ (x , r) = g(T)»i(x) + f(r)i4(x). (6.24) 

Inserting Eqs. ( 16 . 2 1 f) and (I6.22p into Eq. (14. 13f) and using the same parametrization of 
Eq. (16 .4p for the inhomogeneous contributions of the pressure and of the energy density the 
following pair of conditions can be obtained: 

% ll^J + a ^ = ~ aip I 3 + PB ) ' (6 ' 25) 

M = _^l( p (o) +3p (o)). (6.26) 
6 

Supposing, for sake of simplicity, that the barotropic index is constant we shall have that 

„(1) = 1 [PrP + ndrP] 2 PB ((]27] 

P (3w + l) 4a 2 (3w + iy { ' 

where the contribution of the electric fields has been consistently neglected. Inserting now 
Eqs. (I6.2ip and (16. 22ft into Eq. (14.151) the following equations can be readily obtained: 



v {1) - p {1) ■ 

P J Si + Ut(B)-p B 8i 



(6.28) 



d T H + 2H 2 = -^V°> - p (0) ). (6.29) 
Note that, from Eqs. (I6.26P and (I6.29p . the standard Hubble parameter is given by H = 

Hia -3( w +l)/2 where H . 

is an integration constant with the same meaning of the analog con- 
stant introduced in section [51 Using then Eq. (I6.27P to eliminate from Eq. (I6.28p . we 
have that 

tffli i ova ai i d( &r ■ - 1 + w 



dtfi + 2Hd T {3{ + 81 [j^KP + 2 YT3^ ndr(3 
= ^ 2 4(|^)pb - 2a 2 4 - 2aVM£). (6-30) 
Equation (I6.30p can be solved by positing 

P{(x,t) = g(T)f4(x) + /(r)^(x). (6.31) 
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Then, inserting Eq. ( I6.3ip into Eq. ( I6.30P and assuming w 7^ 1/3 we obtain 



9[r) 



v°Ax) 



,3uj+1 



f(r) 
4 



,3«)-l 



H?(3w + 5)(3w + 1) 

4 



J? (£) + 



3w 2 — 6w — 5 



7T 



H?(3w- l)(3w + l) 



BAx)BHx 



4(9w + 5) v 7 

9w 2 - 9w - 8 



9(6u; + 5 — 3u7 2 



-E 2 (f)^ 



(6.32) 
(6.33) 

(6.34) 



where Pl{x) = r- (2, r)a 2 (r) accounts for the intrinsic curvature computed from azij(x) and, 
by definition, 



Bi(x)B j (x) 



a(x) 



B m (x)BHx) 1 B 2 (x) = amn }*> B m (x)B n (x), 

a(x) 



with a(x) = det(ajj). Consequently ■y i j(x,r) can be written as 
If w = 1/3 



a 2 (r) 



aij(x) + ^ij(x)a Zw+1 + Vij(x)c? w 1 



where, in the latter case, 



aij(x) + fiij{x)T 2 + Vij(x) lnr 



3H? 



i4(x) 



f 2 



2nH? 



B 2 (x) 



(6.35) 



(6.36) 



(6.37) 



(6.38) 



If w < —1/3 the contributions of the spatial curvature and of the magnetic field be- 
come progressively subleading even if the terms proportional to are diluted faster. If 
— 1/3 < w < 1/3 the contribution of the spatial curvature is progressively increasing while 
the contribution associated with the magnetic fields decreases. Finally if w > 1/3 both terms 
increase but at a different rate. The solutions illustrated in the present section generalize the 
perturbative treatment where the large-scale magnetic fields are taken as a supplementary 
component in the linearized equations of the relativistic metric perturbations. The linearized 
evolution of the curvature perturbations is usually described not in terms of 7^ but rather 
in terms of an appropriate gauge-invariant combination. It is therefore necessary to consider 
the interplay between the two different descriptions in the case when fully inhomogeneous 
magnetic fields contribute to the curvature perturbations which also receive an independent 
contribution from the standard adiabatic mode. 



7 Nonlinear adiabatic mode in the plasma frame 

The initial conditions or the calculation of the CMB observables in the presence of large-scale 
magnetic fields are set within linear perturbation theory expressed either in a specific gauge 
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or within a suitable gauge-invariant treatment |S]. In contrast with earlier studies, it 
has been recently clarified that different kinds of initial conditions are contemplated ranging 
from the magnetized adiabatic mode to the various magnetized entropy modes. Nonlinear 
generalizations of the magnetized curvature perturbations will now be discussed. 

In the standard perturbative treatment of large-scale inhomogeneities, the curvature 
perturbations on comoving orthogonal hypersurfaces (conventionally denoted by TV) and 
the curvature perturbations on uniform density hypersurfaces (conventionally denoted by Q 
are often used to parametrize the large-scale curvature inhomogeneities. Both variables are 
gauge-invariant and can therefore be expressed in any gauge such as the longitudinal gauge 
or the synchronous gauge (ISJ 0B] . In the longitudinal gauge (see appendix Q3] for an explicit 
definition in terms of the ADM variables) and in the plasma frame TZ and ( can be written 
afi 

K(x,t) = n2 _ m , (7.1) 

mr) = ^ + (7.2) 
3(pt + Pt) 

where 5^p denotes the scalar fluctuation of the energy density of the fluid sources in the 
longitudinal gauge; 6 s pB = -B 2 (i;)/[87ra 4 (r)] denotes the fluctuation of the energy density 
associated with the fully inhomogeneous magnetic field, while p t and p t are, respectively, the 
total energy density and pressure determining the (homogeneous and isotropic) background 
geometry. 

By computing the difference (( — TV) and by comparing the obtained result with the 
Hamiltonian constraint (see Eq. flB.5|) of appendix |B]), the difference (( — TV) becomes 

C~K= \ 7.3 

3 if, a 2 {p t + p t ) 

which is negligible in the limit of typical length-scales much larger than the Hubble radius 
at the corresponding epoch. The variables ( and TZ are gauge-invariant: their numerical 
value does not change for two coordinate systems connected by infinitesimal coordinate 
transformations of the type 

x ^ = + e^, e M = a 2 (r)(e ,-^e) (7.4) 

where both eo and e depend on the space-time pointful. A coordinate transformation of the 
type ( 17. 41) will change metric fluctuations according to the Lie derivative in the direction of 



15 The notation TZ and £ (as opposed to TZ and £) is meant to distinguish the quantities defined in pertur- 
bation theory from their counterpart defined to a given order in the gradient expansion. 

16 As explained in appendix [B] the first part of the present section assumes the standard perturbative 
treatment and, therefore, the underlying background geometry is taken to be conformally fiat as in the case 
of the vanilla ACDM scenario. 
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the vector consider, for sake of concreteness, the transformation from the synchronous 
coordinate system to the longitudinal one: 

^(S) ^ ^(L) = ^(S) + H ^ F (S) ^ F (L) = F (S) _ £) (7 5) 

B (S) B (L) = B (S) +eQ _ (S) ^ (L) = (S) _ _ (7 g) 

In the synchronous gauge < - s - ) = and = while, by definition of longitudinal gauge, 
= and B^ = (see, e.g. [M]). Consequently, the standard relation between the 
variables appearing in ( and 1Z can be written as 

ip = -f + Ud T F, = -d 2 T F - Hd T F, 5 S (L) = 5 S (S) - d T p^d T F. (7.7) 

By inserting Eq. ( I7.7P into Eqs. (JT7TJ) and ( I7.2p . the standard synchronous expression of 1Z 
and ( can be readily obtained 

K = { + ^ w , C = « + f^, (7.8) 
H 2 -d T H 3(p t +pt) 

and it coincides, as expected, with previous discussions (see, e.g. [16]). Because of gauge 
invariance, the evolution equations for ( and 1Z can be derived and discussed in any gauge, 
for instance, by linearizing the covariant conservation of the energy momentum tensor of 
the sources. Let us consider, for sake of concreteness, the matter radiation transition and 
let us define the fluctuation of the total pressure in terms of the sound speed and of the 
non-adiabatic pressure fluctuation 5p na( j: 

&sP = c 2 J s p + 5p nad , 4 = ^ = 4 (7.9) 

a T pt 3 (3a + 4) 

where a = a(r)/a eq and the sound speed, as shown in the case of the matter radiation 
transition, is fully homogeneous. The result for the evolution of ( when the plasma effects 
are carefully included has been given in the second paper of Ref. [4"8] 



a - , nSpmd E-J H5 s p B (3c 2 st - 1) W s p E t (7 m , 

(Pt + Pt) 3a 4 (p t + p t ) 3(p t + pt) 3(p t + Pt) 3 

where the following quantities have been introduced 

^ t) = IW±f« et = v.. t; (,n) 

gi and g 2 are two numerical constants which depend upon the ACDM parameters which can 
be explicitly computed (see second paper of [4"8]). For the fiducial set of ACDM parameters 
g 1 <C g 2 ~ 10~ 5 . Neglecting the electric fields, Eq. (17. lip can be easily solved by direct 
integration since, for length-scales much larger than the Hubble radius, Eq. (17. lip reads: 

gC _ gggnad f 2 1N | 5 spB - „s 

aina (Pt + Pt) V st 3Apt + Pt)' 1 ' J 
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with the result that, in the case of the magnetized adiabatic mode (i.e. 5p na a = 0), the 
solution of Eq. ( 17. 1 2 j) becomes: 

C(x,r) = CM ~ \ R ^)^l- ( 7 - 13 ) 

In the nonlinear case it does not make sense to separate the energy density and the pressure in 
a background value supplemented by the corresponding fluctuations. Still it will be possible 
to define the sound speed in terms of the barotropic index, i.e. 

st[ ' ' d rP + NK( P + P y u } 

Equation (I7.14p coincides with the second relation of Eq. (17. 9p in the fully homogeneous 
case. However, when the description of the geometry is given in fully inhomogeneous terms, 
the two definitions lead to different results. 

Let us now see how the curvature perturbations can be generalized to nonlinear level. A 
possible nonlinear generalization of the the variable ( is given by: 

d = + Vi(p + p B , (7.15) 

3 6G T p 

while a complementary generalization is 

Q = YihVl + Y^±isl. (7 . 16) 

3 Hp + P) 

Note that Eqs. (I7.15P and (I7.16P are equivalent to zeroth order in the gradient expansion 
since d T p = NK(p+p). However, for practical reasons, the definition (I7.16P will be preferred. 
With similar logic, the nonlinear generalization of 1Z can be written as 

V i lnyg V^-Trif 2 ) 
l ~ 3 + 6£UP + P) ' 1 ' 

The variable given in Eq. (I7.15P is inspired by the one defined in Ref. jU] in the case of an 
energy-momentum tensor dominated by a single scalar field. The pair of variables defined 
in Eqs. fl T . 1 6 [) and (17.171) can be compared to their linearized counterpart by expressing the 
extrinsic curvature in a specific gauge. For instance, in the longitudinal gauge, using the 
results of appendix IBI (and, in particular, Eq. ( IB. 2ft ) Eqs. (I7.16P and (17.171) become 



6V>p + 5 sPB 

—yj + 



3(p t +Pt) 



Ki -> V, 



n 2 - o T n 



(7.18) 



which also implies, because of Eqs. (17. ip and (17. 2p . (-^ — > ViC and Tlf^ — > VilZ. By taking 
the difference of Q and IZi and by using the Hamiltonian constraint of Eq. (12. 9p we obtain 
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By expanding the right hand side of Eq. (17.191) to first order in the (longitudinal) metric 
fluctuations, it can be verified that the obtained result coincides with the spatial gradient of 
Eq. (I7.3p . As far as gauge- invariance is concerned, the approach followed here is similar to 
the one invoked in related contexts (see, e.g. [27JII9]): Eqs. (" 17. 16[) and (I7.17P can be shown 
to be gauge invariant to a given order in the gradient expansion. Suppose, indeed, to evaluate 
Eqs. ( I7.16P and (I7.17P not in the longitudinal gauge but in a different gauge, for instance 
the synchronous gauge. If the variables are truly nonlinear gauge-invariant variables they 
must also be gauge-invariant for infinitesimal gauge transformations of the kind discussed in 
Eqs. ( 17.5p -( |?T6l) . With the help of Eqs. ( 1B.9p . the synchronous gauge expression of IZi can 
be written as: 



R (S) 



H 2 - d T % 



iv. 
3 



(7.20) 



(S) 

where the notation 1Z\ reminds that 7Zi is computed in the synchronous coordinate system. 
The first term in squared brackets of Eq. ( I7.20p clearly coincides with the spatial gradient 
of 1Z (see Eq. ( 17. ip ). The second term in squared bracket appearing in Eq. (17.201) breaks 
therefore gauge-invariance but it is of higher order in the gradients. The same kind of 
considerations can also be extended to the case of Q. In summary, recalling Eq. (17.81) (i.e. 
the synchronous form of 7Z and (), we have that, in the synchronous parametrization, 

Cf = V.C + ^^^C^ + ^V 3 ), (7.21) 



nf ] = v^ + -v 



l. 



= n^ + 0(V 3 ). (7.22) 



Equations (17.211) and (I7.22p are rather interesting and suggest that the variables Q and 
IZi defined in Eqs. ( I7.16P and (17.171) can be computed in any suitable gauge. If the two 
gauges differ by an infinitesimal coordinate transformation, the terms leading to a breaking 
of gauge-invriance are of higher order in the gradient expansion. In the specific example of 
Eqs. (I7.2ip and ( I7.22p the extra terms contain three gradients. 

The observation of the previous paragraph fails, however, in the case of more general 
coordinate transformations. Still the results suggested by Eqs. ( I7.2ip and (I7.22p can be 
generalized along the lines of [17J H9] to the case when the coordinate transformations are 
finite but still sufficiently well behaved. Consider, in particular, the following coordinate 
transformation 

X^X'far), T = T(f,r). (7.23) 

As discussed in jHHH], along lines of constant X 1 the old coordinates change as dx^ = d^Tds 
where s is an arbitrary parameter and 

x ' = f ' {:I> + Ia§Wf dT - < 7 - 24 > 
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If the new time coordinate T is non singular the second term in the transformation of Eq. 
(I7.24p can be discarded to leading order in the gradient expansion. 

From Eqs. (I7.16P and (JT7TTJ) the evolution equation of the nonlinear curvature perturba- 
tions can be obtained in explicit terms by using the governing equations and the covariant 
to conservation equation to leading order in the spatial gradients. Let us therefore take the 
first time derivative of Q and let us drop the contribution of the electric fields which are 
subleading to zeroth order in the drift approximation; the result of this simple manipulation 
is: 

- H 9 ^) + suh) v ' (rXf ' +d ^ - ^f 19M] - (7 - 25) 

But to first order in the gradient expansion d T p = NK{p + p) and d T p-s = iNKp^/S. Thus 
Eq. (17.251) becomes: 

Equation (17.261) generalizes Eq. (17.101) as it can be easily appreciated by considering various 
specific limits. Suppose, for instance, that the barotropic index w is constant both in time 
and in space and that pb = 0. Then, c^ t = w and Eq. (17.261) implies d T Q = 0. This is the 
case of the single adiabatic mode in the absence of magnetic fields. 

In general terms the solution of Eq. (17.261) represents a complicated self-consistent prob- 
lem since both c^ t and w will be both space and time dependent; moreover, always from a 
general point of view, Vi(NK) and VjPb can well be of the same order. To simplify the 
situation let us make the following (not completely realistic) assumption 

MJ f < , w{a) = — 7.27 

JS K p B o[ol + 1) 

where, as before, a = a/a eq . Equation (17. 26ft can then be solved directly and the result is: 

CAx, a) — CAx) : r , 7.28 

S V ; SjV ; 4(3a + 4) p v ; 

where the integration constant has been matched with the value of Q determined in the 
linearized approximation. The solution ( 17.281) is interesting but it assumes an expression for 
the evolution of the barotropic index which is only justified in the homogeneous and isotropic 
case. 

A much safer approach for the computation of the curvature perturbations is to first 
obtain the quasi-isotropic MHD solutions at the wanted order in the gradients (as already 
done in section E]) and then to evaluate explicitly IZi and Q. Bearing in mind that, in the 
notations of section El ^7 = a 3 (r)y / a(x)[l + 0(x,r)/2], the explicit expression of 1Zi(x,r) 
can be computed and it is: 

n 



<v r- \ ldia 1 

■CAX, T) = — 



g(r)dip + f(r)diU 



+ 



6(H 2 - d T H) 



d T gdip + d T fdiV 



(7.29) 
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Recalling the explicit expressions of fi(x) and v(x) we have that 



._ . 1 9j« (9jP 3w+1 



4(9w + 7)(9w + l, QO 



3(3w - l)(3w + l)(6w + 5 - 3w 2 )(w + 1) 



^Q B a , (7.30) 



where now fieO?) — ip B 2 (x)/(2AttH 2 ). By using the explicit expression of the Hamiltonian 
constraint the variable Q(x,r) is only sensitive, by construction, to the gradients of the 
magnetic fields, i.e. 



&(x,t) = TZi(x,r) + 



pa 2 (p + p) 

A(Qin 4- 7VQii, 4- 1 ^ 

di£l B a 6w - 1 . (7.31) 



1^ 4(9 W + 7)(9 W + 1) 

6 a 3(3w- l)(3w+ l)(6w + 5-3w 2 )(ti; + 1) 

Both in the case of Eq. (17.301) as well as in the case of Eq. (17.311) the leading order 
result is fully inhomogeneous and it accounts for the large-scale curvature perturbations 
connected, in this specific case, to the adiabatic solution. The spatial gradients affecting 
IZi are induced both by the spatial curvature as well as by the magnetic inhomogeneities. 
Conversely, the spatial gradients of the curvature affect Q only to higher order because of a 
cancellation arising when the solution for TZi is inserted into the fully ingomogeneous form 
of the Hamiltonian constraint of Eq. (17. 19ft . Note that in the case w = 1/3, Eq. (16 .38f) 
implies that v\ is traceless. This implies that the contribution of the magnetic fields to 
the curvature perturbation vanishes to first-order in the gradient expansion in a radiation- 
dominated Universe while the effect of ordinary spatial inhomogeneities does contribute. 
Consequently, it has been shown that while it is difficult to solve explicitly the nonlinear 
generalization of the evolution equation for the magnetized curvature perturbations, the 
same techniques leading to the quasi-isotropic solutions can be used to derive the expression 
of the nonlinear generalization of the magnetized curvature perturbations. 



8 Concluding remarks 

The general relativistic gradient expansion has been combined with the standard tenets of 
the drift approximation used in the description of cold plasmas. Nonlinear effects are typical 
both of general relativistic dynamics as well as of magnetohydrodynamics in flat space-time. 
It is then natural to treat them in a unified perspective where the standard linearized approx- 
imation is not assumed from the very beginning. A fully nonlinear system has been derived 
in terms of the ADM variables and discussed in specific physical limits. The evolution of 
magnetic and curvature inhomogeneities has been treated and solved on the same footing. 
The present findings have also been contrasted with the standard linearized description of 
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magnetized curvature perturbations both in the presence of non-adiabatic pressure fluctua- 
tions as well as in the case of adiabatic initial conditions. After introducing the nonlinear 
magnetized adiabatic mode, an explicit evolution equation for the magnetized curvature per- 
turbations has been derived to leading order in the spatial gradients. While the results of 
the present investigation pave the way for a more thorough scrutiny of nonlinear effects in 
gravitating plasmas prior to photon decoupling, they also fill an existing gap in the present 
literature. Indeed, the treatment of magnetized plasmas usually rests upon equations written 
in homogeneous and isotropic backgrounds supplemented by the corresponding relativistic 
fluctuations of the geometry and of the fluid sources. Here the modest attempt has been 
to suggest and partially develop a description which is independent of the background but 
valid to a given order in the spatial gradients. 
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A Electromagnetic fields and ADM decomposition 

The explicit form of the of the Christoffel symbols in terms of the ADM variables of Eq. 
(12~TD are: 



d T N N k ^V k N N m N n 
r u — _ i 1 k 

N N N " 



L oi - N N 

P 00 = A (3) VVV + d T N { -2NN m K i m + N m{ - z) V m N i 

]\Ti / \ 

- — (d T N + N m ^V m N-N m N n K mn y 

r°. = -—k 

= ( 3 )r* -i K (A 1) 

1 mn 1 mn ' ly mn- K^-^J 

In Eq. (lA.ip the indices are lowered and raised ny using 7^ so that, for instance, K\ = 
"f^Kf-i] ^V; is the covariant derivative defined with respect to the spatial metric 7^ and 
the corresponding Christoffel symbol is: 

i3) ^mn = (-djlmn + d nljm + d mln ^ . (A.2) 

Recalling Eqs. (12. 2p and ( 12. 3 P the traces K, TyK 2 and r are defined as: 

K = Kl TrK 2 = K J i K t j , r = 7 %. (A.3) 

From Eq. (lA.ip the components of the Ricci tensor and the Ricci scalar read 

f) K X7 2 N N m /V« 

= _L_ _ TrA' 2 + _____ _ _ — \7 K + — — £„, (A.4) 

N N N N 

Ri = jjb, (A.5) 
R\ = ±d T Ki-KKi-rl + ±V^N-^V m Ki 

+ ^V m N j K™ - ^ViN m K 3 m - (A.6) 

where, for sake of conciseness, the following notations have been adopted: 

( 3 ) V, = V,, d = (ViK - V k K?) . (A.7) 



32 



From Eqs. (1A.4[) . OA.Sp . flA.6[) and (IA.7[) the components of the Einstein tensors can be 
easily obtained and they are: 
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(A.8) 
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(A.10) 



For a perfect relativistic fluid the components of the energy-moementum tensor with covari- 
ant indices are: 

Too = (p + p)u Q u - p(N 2 - N k N k ), 
T i0 = (p + p)u Ui +pNi, 

T ij = (P + p)UiUj + Pliy (A- 11 ) 



The condition g^ u UnU v = 1 implies, in terms of the ADM variables, 



u - N k u k 



N 2 (l + u 2 ), 



u 2 = j^UiUj. 



(A.12) 



Equations Eqs. flA.110 and (1A.12j) are written in general terms and, therefore, we shall also 
have that, in general: 

-N k u k 



T ° = p+(p + p)u Vl + u 2 
o _ (P + P) 
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T 
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Vl + U 2 Ui, 



NU y/T+ 



Ti = —p5j — (p + pjUiV? — - ^ ' f Vl + u 2 Ui N 3 ' 



u' 
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(A.13) 

(A.14) 
(A.15) 



From Eqs. (1A.9I) and ( 1A.14I) the explicit form of the the Hamiltonian and of the momentum 
constraints, (i.e. Gq = (^Tq and G° = £pT®) becomes, respectively: 



K 2 + r - TiK 2 = 2£ 2 P [p + (p + p)u 2 ), & = t\{p + p)v / TT^ 



Ui. 



(A.16) 



Note that the term proportional to C q in Gq vanishes exactly with the terms containing 
N k u k in Eq. (1A.13I) once the momentum constraint is imposed. In the case Ni = and 
within the parametrization of Eqs. (13.61) and (13.71) the relevant component of the fluid 
energy-momentum tensor can be written as 



rpOO 
rpij 



N 2 



[p cosh 2 y + p sinh 2 y] , 



P + P 
N 



sinh y cosh yv l 



(p + p)v l v 3 sinh 2 y + p^fK 



(A.17) 
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The covariant conservation of the energy-momentum tensor can be written in terms of the 
Christoffel symbols obtained in Eq. ( lA.lj) 



V M T^° = d T 00 + d k T k0 + [2T° 00 + T k ok ) T 00 + [3T 0k + TUT™ + T° kj T kj , (A. 18) 



(A.19) 



V M T^ = d T 0i + d k T ki + T^T 00 + 2Y i 0l T Qj + ( T° 00 + T k k ^j T i0 



+ \- j} ;p k ■ (v% ■ v%yn. 



which can also be explicitly written, in ADM variables, as 

-p + p 



-^{dr[ps 2 (y) + pc 2 (y)] + N 2 d k 



N 



c(y)s(y) v k 



NK[ps 2 (y)+pc 2 (y)} 



3V k N + NTi 



(p + p)s(y)c(y) v k 



- NK 



(p + p)v k Vs 2 (y) +p7^' + n fc3 ' 



(A.20) 



+ 



N 
N 



dr[(p + p)s(y)c(y)v 1 } + d k [(p + p)s 2 (y)v l v k + jrf lk + U 



[pc 2 (y) +ps 2 (y)} 



- 2K)(p + p)s(y) c(y) v j - K(p + p)s(y) c(y) v l 
+ TU(p + p)v k Vs 2 (y)+p^ + U k ^ 



+ 



kjl 
N 



+ T k k )[(p + p)v^s 2 (y)+p^ +n 



(A.21) 



where the shorthand notation c(y) = coshy and s(y) = sinh?/ has been adopted and where 
IP- 7 denotes the possible contribution of the anisotropic stress which has been included for 
completeness. Note that, in Eqs. (1A.20I) and (1A.21I) there is no potential ambiguity since, 
when the shift vector vanishes, T k - coincides with ^T k -. This is not the case in general as 
the last equality of Eq. flA.lj) clearly shows. 

Finally, the components of the electromagnetic energy-momentum tensor become, in 
ADM variables, 



'o 



7-0 

V 

> 1. 



1 

1 



E m E n "v uP • ri ijk 

Ink 1 ijiri'l urn nil 
in D D 

N 2 
QijkE k E j — Pij k B k E j 



Pnm InkV "* ^ \ , -, 
~Jf H 2 ^ ijm ' 



Air^j I V 7p m ^' - la b Pi mn V Jmb B n B a ) + [Q imnlah ^ mh E n B a - p mi — 



(A.22) 

(A.23) 
B m E j 



+ -A 

4 



P m nkla b r] mnb B k B a + ^(p ab B a E b - q ab E a E b 



— QmnklabV mnb E k B a 



(A.24) 
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In Eqs. flA.22p . (IA.23|) and (IA.24|) the following auxiliary tensors have been introduced, 

(A.25) 



Qki 



N k N { + lik (N 2 - N rn N m ) N klijlrnn ^ n 
, Pki = ~ 



v/7 



1 mnp 



/ \ --J .'i. ■ ■ ,jx- d lim r ijn r 1pk'n , * ,-, . 

Ulijk — AT !— , r. ij k — — . [K.ZO) 



In terms of the tensors given in Eqs. ( 1A.25I) and ( 1A.26}) the covariant components of the 
field strengths can be written as: 



Foi — q k iE k — p k iB k , Fij — Q^ k E k — Pij k B k , 



(A.27) 



while the controvariant components of the field strength have been already reported in Eqs. 
(USD and O in terms of £ k = (N/^y)E k and of B k = (N/^f)B k . In the gauge where the 
shift vector vanishes, it is easy to show, within the present decomposition, that 



N ^ 

Qki ~^~)kii 



Pki — 0, Qij k — 0, 
and Eqs. f TA~22l . f TA~23l) and (tOijl reduce to: 



p _ limljnlpk rank 



(A.28) 



'0 

V 



7T7 ' 



47T7 N 



rank 



-B p E j 



47T7 



Si 



lim E m E> + lm B m Bi - -± lmn [E m E n + B m B> 



(A.29) 
(A.30) 
(A.31) 



In terms of the shorthand notation employed in the bulk of the paper, the components 
reported in Eqs. (IA.29I) . ( 1A.30I) and ( 1A.31I) can also be written as: 



r° 

'o 

r° 
v 



i 



57T7 

7i"i 



E-E + B-B\=p B +p E 



47T7 



(E x B) 



Si 



lim E m & + lim B m B> -^(E-E + B-B 



47T7 

-(pv+PB)Si + ni(E) + ni(B), 



(A.32) 
(A.33) 

(A.34) 



where p E = Pe/3 and p& = Pb/3; n{(E) and U-l(B) denote, respectively, the electric and the 
magnetic anisotropic stress. 
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B Relations to Bardeen formalism 



The Bardeen formalism |45] is one the main tools customarily employed for a quantitative 
assessment of the impact of large-scale magnetic fields on the CMB anisotropies. It is 
therefore useful to describe the relation of the methods described in the present paper to the 
Bardeen approach. In the linearized theory of cosmological perturbations the fluctuations 
can be separated in scalar, vector and tensor modes as: 



(B.l) 



The scalar modes of the geometry are parametrized in terms of four independent functions. 
The vector modes are parametrized in terms of two pure vectors Wi and Qi obeying ^\7iQ l = 
and ^ViW % = 0. Finally the tensor modes are parametrized in terms of a rank-two tensor 
hij which is both divergenceless and traceless. Overall, before gauge fixing, the number of 
independent functions amounts to 10. At nonlinear level, the decomposition of Eq. (1B.1I) is 
meaningless but still the more general description discussed in this paper contains, as special 
cases, the gauge dependent approaches to the magnetized CMB anisotropies [8], [9]. For 
instance, the choice of the conformally Newtonian gauge (often dubbed longitudinal gauge) 
corresponds to iV, = 0, N(x, r) = a(r)[l + <p(x, r)] and jij(x, r) = a 2 (r)[l — 2ip(x, r)]5y. The 
extrinsic curvature and ypy are therefore given by 



-- + -[d T ij + H(t) 

a a 



U l 5 



V7 = a 3 (l-3# 



(B.2) 



From Eqs. ( 1A.8j) . ( 1A.9|) and ( 1A.10I) the various components of the einstein tensors can be 
obtained by keeping only the terms which are linear in the metric fluctuations. For instance 
Eq. ( 1A.8|) implies 



G° (X,T 



G^(t) + ^(£,t 



S s G° (x,r) = - 



V 2 ij - 3H(H(j) + ij') 



(B.3) 
(B.4) 



It is useful to remark that, in Eq. ( IB. 41) 
the perturbative approximation. Converse 



.X, T 



< Gq(t) as implied by the validity of 



y, in the gradient expansion, what matters is not 
the absolute magnitude of the perturbation in comparison with the background but rather 
the number of gradients defining the various orders of the expansion. Following the same 
procedure of Eq. (IB. 41) the evolution equations in the longitudinal gauge can be explicitly 
obtained from Eqs. (TA3]) and (jA.lOp : 



1)2 

V 2 ij - miHcj) + d T ij}) = ya 2 [5 s p + S bPb + S sPE ], 



2V 2 (U4> + d T i))+f P 



(Pt + Pt)0t + 



V • (E x B) 

4^ 



(B.5) 
(B.6) 
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dfy + n{d T <f) + 2d T ^) + (H 2 + 2d T H)<j> + -V 2 (0 - V) 



s sP + s aPE + 5 sPB )5f - n? - ni(£) - n|(s) 



(B.7) 



where # t = V ■ #t- In analog terms the evolution equations can be obtained in a different 
gauge either by performing the appropriate gauge transformation on both sides of Eqs. ( IB. 51) . 
(IB.6P and (IB.7[) or by using again the general form of the Ricci (or Einstein) tensors reported 
in appendix [A] The synchronous gauge equations can be obtained by positing, 



N(t) = cl(t), 7y (f , r) = a 2 (r)[(l + 2£)<% + 2<9^F]. (B.8) 
In the case of the synchronous gauge condition the extrinsic curvature and ^Fy read 



m = - 



v / 7 = « 3 (l + 3^ + V 2 F) 



(B.9) 



Finally, not only the evolution of the scalar modes can be readily obtained but also those 
for the vector and for the tensor modes. For instance, from Eq. (1A.6|) . in the gauge 



NAx,t) 



N(t) = a(r), 7^ = a (r)tfy, 
the vector fluctuation of the Ricci tensor can be written 



m(2,T)=l£(T)+6 w MfaT) 



a 2 (r)Qi(x, r) 



(B.10) 



R{( T ) = -*td T H + 2H 2 )5l 



did T Q j + d j d T Q t 



+ 2U 



d .Qi + Q3Q. 



(B.ll) 



In the bulk of the paper we shall not dwell on the derivation of the linearized results. It will 
be however understood that they are easily obtainable from the general expressions reported 
in appendix [A] by going through the same steps outlined in the specific examples sketched 
in this appendix. 
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